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Abstract: The control study of vehicle semi-active suspension with magneto-rheological (MR) dampers has been attracted much 
attention internationally. However, a simple, real time and easy implementing semi-active controller has not been proposed for the MR 
full-vehicle suspension system, and a systematic analysis method has not been established for evaluating the multi-objective suspension 
performances of MR full-vehicle vertical, pitch and roll motions. For this purpose, according to the 7-degree of freedom (DOF) full- 
vehicle dynamic system, a generalized 7-DOF MR and passive full-vehicle dynamic model is set up by employing the modified Bouc- 
wen hysteretic force-velocity (F-v) model of the MR damper. A semi-active controller is synthesized to realize independent control of 
the four MR quarter-vehicle sub-suspension systems in the full-vehicle, which is on the basis of the proposed modified skyhook 
damping scheme of MR quarter-vehicle sub-suspension system. The proposed controller can greatly simplify the controller design 
complexity of MR full-vehicle suspension and has merits of easy implementation in real application, wherein only absolute velocities of 
sprung and unsprung masses with reference to the road surface are required to measure in real time when the vehicle is moving. 
Furthermore, a systematic analysis method is established for evaluating the vertical, pitch and roll motion properties of both MR and 
passive full-vehicle suspensions in a more realistic road excitation manner, in which the harmonic, rounded pulse and real road 
measured random signals with delay time are employed as different road excitations inserted on the front and rear two wheels, by 
considering the distance between front and rear wheels in full-vehicle. The above excitations with different amplitudes are further 
employed as the road excitations inserted on left and right two wheels for evaluating the roll motion property. The multi-objective 
suspension performances of ride comfort and handling safety of the proposed MR full-vehicle suspension are thus thoroughly evaluated 
by comparing with those of the passive full-vehicle suspension. The results show that the proposed controller can ideally improve multi¬ 
objective suspension performances of the ride comfort and handling safety. The proposed harmonic, rounded pulse and real road 
measured random signals with delay time and asymmetric amplitudes are suitable for accurately analyzing the vertical, pitch and roll 
motion properties of MR full-vehicle suspension system in a more realistic road excitation manner. This research has important 
theoretical significance for improving application study on the intelligent MR semi-active suspension. 
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1 Introduction 

Suspension plays an important role in ensuring ride 
comfort and handling safety of the road vehicle and the 
semi-active controllable magneto-rheological (MR) damper 
has been widely explored to realize an intelligent semi¬ 
active vehicle suspension in the past decade. Lots of papers 
were published on modeling the hysteretic force-velocity 
(F-v) characteristic of MR damper and synthesizing semi¬ 
active controller aimed at the MR quarter-vehicle 
suspension 1 ' 7J . Typical models include the Bouc-wen 
hysteretic F-v model proposed by SPENCER, et al [81 and the 
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sigmoid hysteretic F-v model proposed by WANG, et al [9] , 
the former has advantage of precisely describing transient 
process but lacks precision in describing steady state 
property, while the later has high precision in describing 
steady state property but lacks precision in describing 
transient process of the MR damper. The semi-active control 
studies, such as skyhook control, H x control and neural 
fuzzy control, focus on improving the ride performance of 
MR quarter-vehicle suspension by suppressing vertical 
motion vibration. Whereas, the full-vehicle suspension has 
requirements of multi-objective suspension performances 
involving the vertical, pitch and roll motions, few studies 
have considered the semi-active control aimed at MR full- 
vehicle suspension. CHOI, et al [7] , early investigated 
adaptive control of MR full-vehicle suspension with body 
mass variation by applying H x control strategy. YANG, et 
al tl0] , proposed a semi-active controller employing fuzzy 
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T-S neural network strategy for the MR Ml-vehicle 
suspension. LIAO, et al 1 " 1 , proposed a coordinate semi¬ 
active controller combining skyhook, ground-hook and 
fuzzy control schemes. The above achievements have 
strongly put forward the semi-active control study of MR 
Ml-vehicle suspension. However, a proper analysis method 
has not been set up for evaluating the multi-objective 
suspension performances about the vertical, pitch and roll 
motions of Ml-vehicle, and the proposed semi-active 
controllers seem more complicate and difficult in real 
vehicle. 

According to 7-DOF Ml-vehicle dynamic system 1121 , a 
generalized 7-DOF MR and passive Ml-vehicle dynamic 
model is set up by employing the proposed modified 
Boucwen hysteretic F-v model of MR damper, and an 
asynchronous semi-active controller with the proposed 
modified skyhook damping scheme is synthesized to realize 
independent control of the four MR quarter-vehicle sub¬ 
suspension systems in Ml-vehicle. Considering the distance 
between front and rear wheels in Ml-vehicle, the proposed 
harmonic, rounded pulse and real road measured random 
signals with delay time are employed as input excitations on 
the front and rear two wheels from road surface, so as to 
systematically analyze the vertical and pitch motion 
properties, and the above excitations with different 
amplitudes are further employed as the left and right two 
wheel road excitations for evaluating the roll motion 
property, in a more realistic road excitation manner. The 
multi-objective suspension performances of ride comfort 
and handling safety of the proposed MR Ml-vehicle 
suspension are thus thoroughly evaluated by comparing 
with those of the passive Ml-vehicle suspension. 

2 7-DOF Full-vehicle Suspension Model 

It is required to set up a proper MR Ml-vehicle 
suspension model, so as to evaluate the multi-objective 
performances of the vertical, pitch and roll motions. Fig. 1 
shows a 7-DOF MR Ml-vehicle suspension model. 



Fig. 1. 7-DOF dynamic model of MR and passive 
full-vehicle suspensions 


In Fig. 1, Mis sprung mass, Ig and Ig, are inertia moments 
of the pitch and roll motions. x % , 0, 0 represent vertical 
displacement, pitch and roll angles, and a, b and l\, l t are 


distances from the front and rear, as well as the left and right 
wheels to the gravity center of sprung mass, respectively. 
m ui , k si , c si , F MRi , k ti represent unsprung masses, suspension 
stiffness coefficients, suspension damping coefficients, MR 
damper forces, and tire stiffness coefficients, and x si , x ui , x ti 
represent sprung and unsprung masses and road excitation 
displacements of the four quarter-vehicle sub-suspension 
systems, respectively. Herein, the subscript z=l, 2, 3, 4 
denote the four front and rear sub-suspension systems at left 
and right sides of the Ml-vehicle, respectively. The 
proposed suspension model can be easily employed to the 
MR Ml-vehicle suspension by setting c s; =0 and to the 
passive Ml-vehicle suspension by setting F MR ,=0. 

The generalized dynamic equations of MR and passive 
Ml-vehicle suspensions can be thus expressed as follows: 

Mx g = - xj-tkjx si - x ai )-^F MRi , 

l e e = - xj + k si (x si - xj + F MRi - 

b J2 c > j (4 - 4)+4 (** - x ai )+ F MRi , 

V = -AM4 -4)+4(u, 4 0) 

C s3 (4 - 4 ) + 4 <As 3 - *u3 ) + 4 r31 + 

414(4 -4)+40s2 ~ x ui) + f mri + 

C s4 (4 -40+4 (*s4 -*u4) + 4rJ> 

rn m .x ul = c s; (i, -X ui )+k si (x si -x ui )-k tl (x ui -x,) + F mr „ 

where x si and x u; . represent displacement velocities of the 
sprung and unsprung masses, x g , 6 and <j) represent 
vertical displacement, pitch angle and roll angle 
accelerations of the sprung mass, respectively. 

F’mrj represent MR damper forces in the four MR quarter- 
vehicle sub-suspension systems and can be calculated 
through an appropriate hysteretic F-v model of the MR 
damper. In this study, the proposed generalized hysteretic F- 
v model decoupled with current control and hysteresis force 
is employed to modify the proposed original Bouc-wen 
hysteretic F-v model 1131 , consequently, the proposed 
modified Bouc-wen hysteretic F-v model in Eqs. (2)—(4) can 
accurately describe nonlinear hysteretic steady and transient 
characteristics of the MR damper: 

4r,(4av) = C(4)4(l,), (2) 

C(/ d ,) = l +- k? --, (3) 

di ' l+exp (—a 2 (/ di +/ 0 )) l+exp(-a 2 / 0 ) 

F b (v T ) = c l y + k l (x-x 0 ), (4) 

y = — |— [az+c 0 x+k 0 {x-y)\, 
c o+ c i 
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z = -r\x-y\z |z|" _1 -P{x-y)\z\' + A(x - y), 

where C(z d ,) formulates the MR damper force in relationship 
with driving direct current (z' di ), which can accurately 
describe the nonlinear variation and saturation properties of 
magnetic field depending on the controlled driving current. 
F h (y„) formulates the MR damper force in relationship with 
relative motion velocity ( v ri = Jc si — x aj ) of damper piston, 
which is expressed with the Bouc-wen hysteron. 

A candidate MR damper with working source 12 V and 
limited maximum drive current 0.5 A, manufactured by 
CARRERA company 13-51 , is employed in this study, and the 
tested data are used to identify parameters of the proposed 
modified Bouc-wen hysteretic F-v model, such as 
^=8 816.9, p=233 849.1, a=20 373.7, c 0 =l 368.7, Cl = 
6 222.7, Ab=184.1, k x = 1 528.1, 4=10.092, n= 2, ,4=20.6, 
Xq= —0.004, a 2 =7.526, 7 0 =0.069. Fig. 2 shows comparison 
of the model results and measured data under varied drive 
current (0-0.4 A) and harmonic excitations with amplitude 
12.5 mm at frequency 1.5 Hz and 2.5 Hz, respectively. The 
result illustrates that the proposed modified Bouc-wen 
hysteretic F-v model can precisely describe real operation 
characteristic of the MR damper 1131 . 



-1 500l_i_i_i_i 

-0.10 -0.05 0 0.05 0.10 

Relative motion velocity v r /(m-s“') 

(a) / d =0—0.4 A, a m =12.5 mm, /= 1.5 Hz 



(b) i = 0-0.4 A, a m =12.5 mm, f- 2.5 Hz 


Fig. 2. Result comparison of model calculation 
and measured data 


3 Skyhook-based Semi-active Controller 
of MR Full-vehicle Suspension 

The original skyhook damping scheme, proposed by 
KARNOPP, et al 1141 , intends to suppress vibration 
acceleration of sprung mass by increasing and decreasing 
the suspension damping coefficients when the sprung mass 
velocity is in same and inverse phases with relative velocity 


of the damper piston, respectively, so as to ideally improve 
the ride comfort suspension performance. A modified 
skyhook semi-active control strategy in current form for the 
MR quarter-vehicle suspension, proposed by the authors 13-51 , 
is directly applied to synthesize an asynchronous semi¬ 
active controller in Eq. (5), which is helpful to greatly 
simplify controller design complexity of the MR full- vehicle 
suspension, and can realize independent control of the four 
MR quarter-vehicle sub-suspension systems. The proposed 
modified skyhook-based asynchronous semi-active 
controller has further merits of easy implementation in real 
application, wherein only absolute velocities with reference 
to the road surface of sprung and unsprung masses are 
required to measure in real time when the vehicle is moving. 

4 = kk,l', ,. = 1>2i3>4> (5) 

I 0 , 4 . (jr s ,.-i ul . )<0, 

where z d/ are controlled drive currents of the four MR 
dampers, x si and x ui are displacement velocities of sprung 
and unsprung masses in the four MR quarter-vehicle sub¬ 
suspension systems, respectively. 4 denotes controller gain 
and m(m> 2) denotes controller order, which has potential 
to suppress the hysteresis effect and to improve the 
damping force bandwidth. In this study, the controller 
parameters are chosen as m=2 and 7 d =3, and the 
parameters of 7-DOF full-vehicle model of the Ford 
Granada car are listed in Table l 112 l 

Table 1. Parameters of the Ford Granada ear 


Parameter Value 

Sprung mass of vehicle M/kg 1 480 

Inertia of pitching motion 7 d /(kg • m 2 ) 2 440 

Inertia of rolling motion /®/(kg ■ m 2 ) 380 

Front unsprung mass m„\, m^i/kg 40.5 

Rear unsprung mass m^m^/kg 45.4 

Tire stiffness coefficient k t i-k, 4 /(kN • m" 1 ) 192 

Front suspension stiffness coefficient k s i, k s2 /(kN • m -1 ) 17 

Rear suspension stiffness coefficient k^k ^/(kN • m" 1 ) 22 

Suspension damping coefficient Cji-c^/CkN ■ s • m _1 ) 1 500 

Distance from front suspension to center of gravity a/m 1.25 

Distance from rear suspension to center of gravity b/ m 1.51 

Half of vehicle width 4, /,/m 0.74 


4 Analysis of Multi-objective Suspension 
Performances 

The proposed skyhook-based asynchronous semi-active 
controllers are applied to the four MR dampers of 7-DoF 
Ford Granada dynamic system shown in Fig. 1, and a 
simulation platform with Matlab/Simulink is built up for 
analyzing both MR and passive full-vehicle suspension 
performances of the Granada car 1151 . Considering delay time 
property between the front and rear two wheels when the car 
is moving, the proposed harmonic, rounded pulse and real 
road measured random signals 13-51 are modified with delay 
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time (t 0 ) to imitate different road excitations, so as to 
analyze steady and transient state properties as well as 
multi-objective suspension performances of vertical, pitch 
and roll motions of the MR full-vehicle suspension in a 
systematic analysis manner. Herein, moving speed of the car 
is assumed as v=60 km/h, and delay time of the excitations 
inserted on the two rear wheels than those inserted on the 
two front wheels can be thus calculated as 
#o=(a+&)/v=0.17 s. 

4.1 Harmonic excitations 

The harmonic excitation is usually applied to analyze the 
steady state response quality of a system. The ride comfort 
suspension performance can be evaluated by analyzing the 
steady state time responses of vehicle body gravity center 
displacement acceleration ( a g = x g ) and unsprung mass 
displacement accelerations ( a ui = x ui ) of the four quarter- 
vehicle sub-suspension systems under harmonic excitation 
at chosen frequency around natural frequency of the sprung 
mass, and by analyzing the frequency responses of 
displacement acceleration transmissibility (Tag) of vehicle 
body gravity center and dynamic load coefficients of the 
four wheels (DLC,) under varied amplitude harmonic 
excitations, in which Tag and DLC t are defined as ratios of 
Root Mean Square (RMS) of a g to RMS of the excitation 
acceleration, and RMS of each wheel dynamic force 
( F u =k u {x ai —Xii) ) to the wheel static force under a given 
harmonic excitation, respectively 13 5 l Moreover, the 
handling safety suspension performance can be reasonably 
evaluated through analyzing the steady state time responses 
of pitch and roll angular accelerations of vehicle body 
motion, under asymmetric amplitude harmonic excitations 
inserted on the left and right wheels. 

The vertical motion property is analyzed by employing 
the harmonic excitations with delay time. Fig. 3 shows the 
steady state response comparison of vehicle body mass 
gravity center displacement acceleration a g and unsprung 
mass displacement acceleration a ul of front sub-suspension 
at left side in both MR and passive full-vehicle suspensions, 
under a harmonic excitation with amplitude 2.0 cm at 
frequency 1.5 Hz, in which the excitation inserted on the 
rear two wheels is delayed with 0.17 s than that inserted on 
the two front wheels. The results illustrate that the peak 
amplitude of a g in the MR full-vehicle suspension decreases 
30% than that in the passive full-vehicle suspension, which 
yields potential effectiveness for suppressing resonance of 
the sprung mass, whereas the peak amplitude of a u i in the 
MR full-vehicle suspension increases 100% than that in the 
passive full-vehicle suspension, which sacrifices part of 
resonance suppression suspension performance of the 
unsprung mass. Moreover, a ul yields obvious self-excited 
oscillations due to the strong hysteresis and force saturation 
nonlinear properties of the MR damper, and yields 
undesirable abrupt responses at about 1.05 s, 1.40 s, 1.75 s, 
due to discontinuous logic condition of the semi-active 
control manner as shown in Eq. (5). 




Fig. 3. Steady state response comparison of MR and passive 
full-vehicle suspensions under harmomc excitations 
with delay time 


Furthermore, the pitch and roll motion properties are 
analyzed by employing the harmonic excitations with delay 
time and asymmetric amplitudes. Fig. 4 shows the time 
response comparison of pitch angular acceleration ( a e = 6 ) 
and roll angular acceleration ( = <j> ) of vehicle body 

motion, herein, in order to conveniently evaluate the roll 
motion suspension performance, asymmetric harmonic 
excitations with different amplitudes of 2.5 cm and 0.5 cm 
at same frequency 1.5 Hz are inserted on the left and right 
two wheels, respectively. Meanwhile, the excitations 
inserted on the two rear wheels are delayed with 0.17 s than 
those inserted on the two front wheels. The results show that 
amplitudes of a 0 and a® in the MR full-vehicle suspension 
are obviously lower than those in the passive full-vehicle 
suspension, which illustrate that the proposed skyhook- 
based semi-active MR suspension has potential benefit of 
controlling the motion posture of vehicle body, and thus 
enhances the handling safety suspension performance. 
Furthermore, the results also show undesirable self-excited 
oscillation and abrupt responses, due to hysteresis of the 
MR damper and discontinuity of the semi-active control 
manner. 

Fig. 5 shows the frequency response comparison of 
displacement acceleration transmissibility T ag of vehicle 
body mass gravity center and dynamic load coefficient 
DLC\ of front wheel at left side in both MR and passive 
full-vehicle suspensions, herein, the harmonic excitation 
with varied amplitudes as shown in Eq. (6) is employed to 
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yield constant amplitudes at lower frequencies (/ < f T ) 
and constant accelerations at higher frequencies (/ > / T ), 
so as to avoid operation saturation of MR suspension 
dynamic system at higher frequencies due to nonlinear 
property of the MR damped 3 51 . 


> m sin(27t/i), f<ft; , = 1234 

aJJ 7 /f)sm(2nft), f>f T , 


( 6 ) 


The above function parameters are chosen as the amplitude 
a m =1.0 cm, the cut-off frequency / T =l .3 Hz, and the 
frequency/=0-25 Hz. 


-Passive 



-Passive 



(b) Roll angular acceleration 

Fig. 4. Steady state response comparison of MR and passive 
full-vehicle suspensions under harmonic excitations 
with delay time and asymmetric amplitudes 



(a) Displacement acceleration transmissibility 


of vehicle body gravity center 



(b) Wheel dynamic load coefficient 
Fig. 5. Frequency response comparison of MR and passive 
full-vehicle suspensions under harmonic excitations 
with varied amplitudes 


4.2 Rounded pulse excitations 

The moving vehicle often encounters harsh road surfaces 
such as hills and hollows, which can be formulated as the 
rounded pulse excitation to evaluate the shock attenuation 
suspension performance of vehicle suspension 13-51 . Herein, 
the rounded pulse excitations inserted on the two front 
wheels are generated using Eq. (7), and the excitations 
inserted on the two rear wheels are considered with a delay 
time in Eq. (8): 


x a, 2 = 0.25a m exp 2 {/uco 0 t) 2 exp (—jua> 0 t), (7) 


The results shown in Fig. 5 indicate that the peak value of 
Tag around 1.5 Hz in the MR suspension decreases 40% than 
that in the passive suspension, while Tag yields slight 
increases in the middle frequency range (2.0-8.0 Hz), which 
further illustrate that the proposed MR full-vehicle 
suspension has powerful role in suppressing the sprung 
mass resonance but sacrifices part of vibration isolation 
suspension performance. The results further show DLC\ 
yields slight increases in the middle frequency range (2.0- 
6.5 Hz) and slight decreases in the higher frequency range 
(6.5-15.5 Hz), which is helpful to improve the road holding 
suspension performance of vehicle, and the sprung mass 
resonance frequency of MR suspension yields a little shift 
than that of the passive suspension. 


x ;3 4 = 0.25a m exp 2 (jua> 0 (t - 1 0 )) 2 exp (-jua> Q (t - t Q ))u{t - 1 0 ), 

(8) 

where the function parameters are chosen as amplitude 
a m =2.0 cm for the left two wheel excitations and a m = 1.0 
cm for the right two wheel excitations, fundamental 
harmonic frequency fflb=10.4 rad/s, pulse stiffness ju=3, 
and delay time f 0 =0.17 s. 

Fig. 6 shows the transient state response comparison of 
vehicle body gravity center displacement acceleration a s , 
pitch angular acceleration a g , roll angular acceleration a 0 
and tire dynamic force F tl of front wheel at left side, herein, 
the wheel dynamic force is defined as the dynamic load of 
each wheel, expressed by F ti =k u {x ui —Xii). The results 
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generally show that the transient processes of a g , a 0 , a,„ and amplitudes of a%, a e and a 0 in the MR suspension yield a 
F tl in MR suspension are greatly shorten than those in the little increases than those in the passive suspension, which 
passive suspension, which benefits to improve the shock sacrifices part of ride comfort and handling safety 
attenuation suspension performance. Whereas, the peak suspension performances. 






Fig. 6. Transient process response comparison of MR and passive full-vehicle suspensions under rounded pulse excitations 
with delay time and asymmetric amplitudes 


4.3 Random excitations 

The real excitation inserted on wheels is usually random 
signal arising from road surface roughness when vehicle is 
moving on a road, and the real road measured roughness 
data reported in Canada 13-5,161 are employed in this study, so 
as to actually evaluate the proposed skyhook-based 
asynchronous semi-active controller of MR full-vehicle 
suspension. The employed real road measured roughness 
data was processed by assuming that the moving speed of a 
car is v=50 km/h, such that the delay time of the excitations 
between the rear and front two wheels is to=(a+b) /v=0.199 
s. Furthermore, amplitudes of the excitations inserted on the 
two right wheels are decreased 50% than those inserted on 
the two left wheels, so as to comprehensively evaluate the 
pitch, vertical and roll motion properties. Herein, the power 
spectral density(PSD) responses 1161 of a g , a 0 , a 0 and F tl are 
analyzed, and the results, as shown in Fig. 7, show that the 
magnitudes of PSD of a g , a s and F a in MR suspension 
greatly decrease in the low frequency range(0.5-2.0 Hz) 
than those in the passive suspension, and thus illustrate that 


the proposed skyhook-based asynchronous semi-active 
controller has potential of the resonance attenuation of 
sprung mass, whereas a$ shows little variations. The results 
further show slight increases of the magnitudes of PSD of 
u g , ag , a0 and F t \ in MR suspension than those in the 
passive suspension in the middle frequency range (2.0-8.0 
Hz), which illustrate that the improvement of resonance 
attenuation suspension performance of sprung mass will 
sacrifice part of the resonance attenuation suspension 
performance of unsprung mass. 

5 Conclusions 

(1) According the modified Bouc-wen hysteretic F-v 
model of the MR damper proposed by authors, the 
established 7-DOF dynamic model comprising MR and 
passive full-vehicle suspensions puts forward an effective 
platform for the semi-active control study of MR full- 
vehicle suspension. 

(2) Considering the distance between front and rear wheels 
in the real vehicle, the proposed harmonic, rounded pulse and 
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real road measured random excitations with delay time and 
asymmetric amplitudes set up a systematic analysis method 
for evaluating the multi-objective suspension performances of 
vertical, pitch and roll motions of MR full-vehicle suspension. 






Fig. 7. PSD response comparison of MR and passive 
full-vehicle suspensions under real road measured random 
excitations with delay time and asymmetric amplitudes 


(3) Employing the proposed modified skyhook-based 
semi-active control strategy aimed at MR quarter-vehicle 
suspension, a skyhook-based asynchronous semi-active 
controller of MR full-vehicle suspension is synthesized for 
realizing independent control of the four MR quarter-vehicle 
sub-suspension systems, which has potential to 
comprehensively improve multi-objective suspension 
performances of the ride comfort and handling safety of 
full-vehicle. 

(4) The above achievement establishes a solid fundamental 
for further study on decoupling semi-active control of the 
four quarter-vehicle suspensions in the MR full-vehicle 
suspension system. 
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